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ABSTRACT 

This paper presents the results of a study of the Barnes 
static Earth sensor assembly (ESA) using on-orbit data 
collected from the Tropical Rainfall Measuring Mission 
(TRMM) spacecraft. It is shown that there exist strong 
correlations between the large penetration angle residuals 
and the voltages produced by the Offset Radiation Source 
(ORS). It is conjectured that at certain times in the 
TRMM orbit the ORS is operating out of its calibrated 
range, and consequently corrupts the penetration angle 
information observed and processed by the ESA. The 
observed yaw drift between Digital Sun Sensor (DSS) 
observations is shown to be consistent with predictions by 
a simple roll-yaw coupling computation. This would 
explain the large drifts seen on TRMM, where the 
propagation of the yaw angle between DSS updates does 
not take into account the possibility' of a non-zero roll 
angle error. Finally, the accuracy of the onboard 
algorithm used when only three of the four quadrants 
supply valid penetration angles is assessed. In terms of 
procedures used to perform this study, the analysis of 
ESA penetration angle residuals is discovered to be a verv 
useful and insightful tool for assessing the health and 
functionality of the ESA. 

INTRODUCTION 

This study was initiated in direct response to several types 
of anomalous behavior observed on the Tropical Rainfall 
Measuring Mission (TRMM) since its launch on 
Thanksgiving Day, 1997. The anomalous behavior 
includes: (1) occasional pixel jumps observed by TRMM 
scientists; 1 * (2) a region of every orbit in which large 
penetration angle residuals are observed; (3) occasional 
occurrences of very large pitch and roll pointing errors 
during periods of Sun interference in one of the quadrants 
of the Barnes static Earth sensor assembly (ESA); and (4) 
larger than desired yaw drifts between yaw updates from 
the two Digital Sun Sensors (DSSs). Although the 
spacecraft is still controlled to within the pitch and roll 
attitude accuracy specification, large pitch and roll 
pointing errors affect the integrity of one of the science 
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instruments and, as will be shown in this paper, may be 
responsible for large drifts in the yaw angle, which 
occasionally exceed the 3cr yaw accuracy requirement. 

The TRMM Barnes static ESA is composed of four 
relatively independent quadrants measuring attitude errors 
about one of the ESA axes.— 1 Therefore, opposite 
quadrants form a redundant pair of attitude sensors, and in 
principle (as was recently done for AQUA J ) it is 
sufficient to process measurements from only one pair of 
adjacent quadrants. The TRMM attitude control system 5 
takes advantage of this redundancy by usina 
measurements from only three quadrants during quadrant 
interference by the Sun or the Moon. One of the purposes 
of this paper is to study the accuracy of this three- 
quadrant algorithm in a more systematic fashion bv 
comparing the resultant roll and pitch pointing errors to 
those obtained using all four quadrants. The accuracy of 
the onboard three-quadrant algorithm is also compared 
w ith that of a more general algorithm, which can handle 
the case when only two adjacent quadrants are available. 

The paper utilizes a new approach, which analyzes ESA 
penetration angle residuals. This analysis is found to be a 
very* useful and insightful tool for assessing the health and 
functionality of the ESA. 

MISSION OUTLINE 

TRMM is an Earth science mission designed to measure 
the rate of rainfall and the total rainfall occurring over the 
tropics and subtropics (between the north and south 
latitudes of 35 degrees). TRMM is a three-axis stabilized, 
Earth-pointing spacecraft, with the spacecraft’s body z B 
axis directed nominally toward the geodetic nadir. The 
Earth pointing (rollpitch) is controlled by the attitude 
control system (ACS) based on measurements taken by a 
Barnes Earth sensor assembly.- 3 6 - 10 (To be accurate, this 
sensor ^ measures deviations from a "horizon bisector” 
vector." 6 The difference between the horizon bisector and 
geodetic nadir vectors is ignored in this paper since it has 
been shown 1 6 that this difference is less than 0.01 degree 
for TRMM.) The yaw angle is updated twice per orbit 
when the Sun vector crosses the x 3 y 8 body plane in the 
field of view (FOV) of one of the DSSs. It is propagated 
between the two yaw updates using yaw gyro 
measurements under the assumption that the ~ 8 body axis 
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Figure 1. Schematic Configuration of the Static Earth Sensor Assembly 

Fig. 1.) The A and B detectors directly measure the 
angular penetration of the Earth limb in their FOV, while 
the S detector views cold space and is used as a reference. 
Three raw voltages — V^j (a single signal from both Aj 
detectors), V Bo , and V S j — measured by detectors A J? 
Bj, and Sj (j= 1, 2, 3, or 4) are converted by the 
spacecraft onboard computer (OBC) into a ‘penetration 
angle’, X } , which is measured from the base of the B } 
detector to the horizon point 

Xj-Pj-Sj, (1) 

where the radial mounting angle 5j is measured from the 
Ze axis to the base of the detector B r All the 5^ 
nominally have the same value, 5 , which is a function of 
the nominal mission attitude (5 = 70.4 degrees for 
TRMM). 


is pointed precisely along the nadir vector. 

THE BARNES STATIC ESA 
Operational aspects of the ESA are shown schematically 
in Fig. 1, with the ellipse denoting the Earth infrared 
horizon. It is composed of four detector assemblies 
referred to below as quadrants Q { - Q 4 . It is convenient to 
introduce 3,6 *' 10 the nominal ESA coordinate system, x £ y E 
ze, in such a way that opposite quadrants lie in either the 
x £ ze or y £ z £ plane. Let x 3 y B z a denote the TRMM body 
coordinate system (BCS) with z 3 nominally pointing 
toward the nadir. The BCS and nominal ESA coordinate 
systems are related through a rotation <t> = 135 degrees 
about the common (z B - z £ ) axis. The quadrant observes 
the horizon point ? } with co-elevation pj measured from 
the Ze axis. It is comprised of four infrared detectors, 
each with a triangular field of view (FOV) of height 5.2 
degrees, two labeled A Jf and one each labeled Bj and S r 
(For clarity, detectors are labeled only in quadrant Q 3 of 
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Most of the voltage from each detector comes from the 
radiation emitted by the detectors to space, so that 
changes in the horizon image appear as small changes in 
large voltage outputs. To detect these changes, the outputs 
from the S-detectors are used to subtract out the signals 
generated in the horizon viewing detectors. If all 
detectors were perfectly matched, no other corrections 
would be necessary. Realistically the match is not 
perfect, and the so-called Offset Radiation Source (ORS) 
is used to compensate for heat lost to space. The source 
replaces the radiation emitted by the detectors to space 
and thereby reduces the large offset signals. The radiation 
received from the ORS by each of the detectors should be 
approximately the same. It has been recognized by 
Barnes Engineering Company that for different 
manufacturing reasons the detector outputs cannot be 
matched to the required accuracy and some additional 
optical adjustment means are provided. An analysis of 
TRMM flight data performed below indicates that this 
optical adjustment is not yet sufficient. As a result, 
compensation for heat lost to space slightly differs for 
different quadrants, which results in noticeable pointing 
errors. 

Reference values of penetration angles can be 
conveniently represented as follows 6 ’ 9,10 


the spacecraft latitude, respectively, and the quadrant- 
T 

dependent term, Ap j , has the form: 


T l 

Ap J --T asin ( 2 Pc) 


nIx, 

j = l,2 

N E.y> 

j = 3,4 


( 6 ) 


where N E>X and N E , y are the projections of the North Pole 
unit vector N on the axes x E and y E . 

Under nominal conditions the TRMM velocity vector is 
either parallel or anti-parallel to the body axis x 3 (the so- 
called ^-forward' and ‘-x-forward’ modes) so that the 
oblateness correction (5) takes the form 



Ap D - k/2 , 

j=l,2 

(7a) 


Ap D +K/2, 

i — , . 

II 

(7b) 


where 

Ap D = --!-asin(2p c )(tan 2 p c sin 2 X + cos~ X) , (8) 

4 

(cf. Eq. (29) in Ref. 11) and 9,10 


v ref 

x j 


x n°m _ (_j)j j = 12 

Xj nom -(-l)A j = 3,4 


( 2 ) 


where 0 and <|> are Earth pointing errors in the ESA 
coordinate system, 


nom 

J 



( 3 ) 


and p® is the nominal angular separation between the 

geodetic nadir and the boresight of the j* quadrant. As 
initially shown by the author, 9,10 the oblateness correction 
to the Earth infrared radius p c , 


% D _ D „ 

* Pj - Pc 


(4) 


— a cos / cosC 
2 s 


sin(2p c ), 


(9) 


where / is the orbit inclination, and ^is the angle between 
the spacecraft velocity and the North Pole unit vector N . 

For TRMM, the quadrant-dependent correction (9) does 
not exceed 0.06 degree and can be neglected when 
compared with observation errors; 12 * however, for the 
sake of generality we will keep this term in the following 
formulas. 


When all four quadrants work, Earth pointing errors, 9 
and <f>, in the ESA coordinate system can be computed 
using the differences of penetration angles measured by 
opposite quadrants 2 ' 5 


can be represented as 


0 * (X,-X : )/2, (10a) 

<j> * (X 3 - X 4 ) /2 _ (10b) 


Ap^ * — - a sin 2 X sin(2p c ) tan 2 p c + Ap ' (5) 

J 4 J 

where a and X are the Earth flatness coefficient 6,1 1 and 


This explains the choice of the label 9 or (J> for the 
quadrants Q { - Q 4 in Fig.l. The symbol *+* or near 9 
in Fig. I indicates whether the measured horizon point is 
located on the positive or negative side of the x £ axis, 
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and similarly for <j) and the y E axis. 

Tn principle, any pair of adjacent quadrants can be used to 
compute 9 and <j>. If our knowledge of the Earth infrared 
radius were precise, each pair of Eqs. (2) would give the 
same answer as Eqs. (10a) and (10b). Due to errors of 
various kind, this is not true. One may expect that 
processing errors would be at least partially eliminated 
if the pointing errors, 9 and <j>, are computed using 
Eqs. (10a) and (10b), which do not require knowledge of 
the Earth infrared radius. The latter formulas, however, 
become inapplicable in the case of quadrant failure, so 
that one has to use the original formula (2) for the pair 
containing the failed quadrant. Below we refer to this 
scheme as the general method (GM) because it applies 
even if just two adjacent quadrants are working. (Note 
that the GM always uses Eq. (10a) or (10b) to compute 
the appropriate pointing error 9 or <j>, respectively, for the 
pair containing two working opposite quadrants.) 

In the case of the failure of a single quadrant, there is an 
interesting alternative, 2 - 3 which explicitly exploits the fact 
that the body x B and y s axes are rotated over the angle of 
135 degrees relative to the ESA x E and y,. : axes. The roll 
and pitch pointing errors, r and p, are proportional to the 
sum of <j> and 9 and to their difference, respectively, 

r ~ + ©) , (1 la) 

p ~^(<J>- 0 )- (lib) 

In the three-quadrant case there are always two pairs of 
working adjacent quadrants. One of these pairs can thus 
be used to express $ + 9 in terms of the difference 
between two penetration angles, whereas <j> - 9 can be 
expressed in terms of the difference between penetration 
angles measured by the second pair of adjacent quadrants. 
Therefore (depending on which quadrant has failed), roll 
and pitch pointing errors can be computed by selecting 
the appropriate formula from each pair of equations: 

r *7f(X:-X j+k). (12) 

V ~ 

or 

r = -^(x 4 -X,-K) (12') 

and 

P 55 ^ C^3 “ ^ I ~ K ) » (13) 


or 

p.-jL(X 2 -X 4+ K). (13') 

Eqs. (12), (12'), (13), and (13') are referred to below as 
the adjacent quadrant method (AQM). It is used by the 
TRMM OBC control system 4 5 to estimate roll and pitch 
angles during interference of a single quadrant by the Sun 
or the Moon. One major advantage of the AQM over the 
GM is that it does not require knowledge of the Earth 
infrared radius, and is therefore expected to be more 
accurate than the direct formulas (2) used by the GM. 

ANALYSIS OF TRMM FLIGHT DATA 

This section describes the results of an analysis performed 
using the TRMM Ground Attitude Support utilities. 8 ’ 1 ' 1 
The TRMM attitude was estimated using measurements 
from two DSSs, ESA penetration angles, and gyro data 
and solving simultaneously for the epoch attitude and 
gyro biases. Since roll, pitch, and yaw angles are 
commanded to zero by the TRMM control system, their 
estimated values are referred to below as pointing errors. 
In most cases two parallel runs were made using “scalar” 
and ''vector” models for ESA measurements. The loss 
function of the “scalar” measurement model is 
constructed using penetration angle residuals for all four 
quadrants. The “vector” measurement model minimizes 
the magnitude of the nadir vector residual, and the 
projections of the vector residuals on the x a and -y 3 axes 
are referred to below as pitch and roll measurement 
errors. (Since the OBC nullifies the observed projections 
of the nadir vector on these axes, the aforementioned 
projections of the vector residuals differ only by sign 
from the corresponding projections of the reference nadir 
vector which explains the choice of the axis directions in 
the definition of the pitch and roll measurement errors.) 

Figs. 2a and 2b depict nominal and measured values of 
penetration angles (X™ m and X™ 35 , respectively) in 

the +x-forward nominal mode. Note that each pair of 
opposite quadrants is represented by a single curve in 
both figures, but for different reasons. The nominal 
values of penetration angles depicted in Fig. 2a were 
computed under the assumption that there were no Earth 
pointing errors, whereas the penetration angles measured 
by opposite quadrants coincide simply because the 
TRMM ACS commands the spacecraft attitude to 
eliminate any differences in penetration angles for each 
pair of opposite quadrants. 

A striking feature of Fig. 2b is a deep minimum near the 
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Figure 2. Anomalous Behavior of Measured Penetration Angles Near the Orbit Apogee 


orbit apogee (see the altitude plot in Fig. 2c). As one may 
notice from Fig. 2a, nominal penetration angles also have 
smaller values near the apogee (due to a slight decrease in 
the Earth angular radius p c ), but the minimum in Fig. 2b 
is significantly deeper. Anomalously small values of 
penetration angles measured by the quadrant Q 3 are 
accompanied by noticeable errors in roil and pitch 
pointing, which implies that the anomalous behavior of 
penetration angles cannot be attributed to inaccuracies in 
the Earth modeling. (Remember that both the onboard 
control and the “vector” measurement model [used to 
generate ground attitude solutions] utilize Eqs. (10a) 
and (10b) and therefore are unaffected by any errors in the 
Earth modeling.) We thus conclude that the observed 
anomalous behavior of penetration angles are caused 
during the processeing of measurement errors. 

The most surprising feature of this anomalous pattern is 
that penetration angle residuals in the quadrants Q3 and Q 4 
increase simultaneously (see Fig. 3a). One would expect 
that all local changes in the Earth infrared angular radius 
(such as horizon radiance variations 1 : ) should affect only 
one of the quadrants. There must be a coupling 


mechanism that simultaneously affects voltage outputs 
from all four quadrants. The ORS, an onboard heat 
source that simultaneously compensates all four quadrants 
for heat lost to space, thus becomes a candidate for the 
source of this anomaly. 

Fig. 3 reveals a pronounced correlation between large 
penetration angle residuals and ORS voltages. As seen 
from the figure, the anomalous behavior of penetration 
angles happens near the maximum of the ORS 
temperature, which results in a noticeable decrease in the 
ORS voltage. It is unclear, however, whether the observed 
errors are directly caused by the increase in the ORS 
temperature. In fact, both the decrease in the ORS 
voltage and the accumulation of errors in penetration 
angles reverse themselves after the spike, which happens 
7 minutes before the ORS temperature reaches its 
maximum. This behavior remains unexplained. 

In the example discussed above the anomalous behavior 
of penetration angles happen to occur near the orbit 
apogee, but it equally often occurs near the perigee, as 
illustrated by Fig. 4, which also uses data for the +x- 
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Figure 4. Anomalous Behavior of Measured Penetration Angies Near the Orbit Perigee 


6 


American Institute of Aeronautics and Astronautics 









AIAA-2000-4245 



Time in seconds since 19990722.000000 


Figure 5. Penetration Angle Residuals for ‘Regular’ ORS Temperatures 


forward mode. The curve representing penetration angles 
measured by the third quadrant now has the second 
minimum, which occurs near the apogee and which is 
again accompanied by relatively large Earth pointing 
errors. 

As seen in Fig. 3, penetration angle residuals start 
growing when the ORS temperature exceeds 20. 5C°. For 
this reason ORS temperatures lying within 20±0.5 C° will 
be referred below as ‘regular’. Fig. 5 magnifies the 
fragment of Fig. 3a associated with ‘regular’ ORS 
temperatures. Note that penetration angle residuals for 
the quadrants Qi and Q : have opposite sign. As a result, 
use of three-quadrant formulas (12) and (13') would 
reduce roll and pitch errors, compared with the 
conventional four-quadrant algorithm. This is obviously a 
simple coincidence, but it shows that the Earth modeling 
is done sufficiently accurate and there is no clear 
preference between different formulas. In particular, 
sufficiently accurate estimation of roll and pitch pointing 
errors can be done using measurements only from two 
adjacent quadrants. 


life examples of quadrant interference, when the TRMM 
control system has to use the AQM for Earth pointing. 
Fig. 6 depicts a period in which QI experiences Sun 
interference. The interference period is marked by two 
vertical dotted lines in Fig. 6b. There is practically no 
change in the penetration angle residuals (and therefore in 
pointing errors) immediately after the interference period 
begins, and where the TRMM onboard control system has 
stopped processing voltages from the first quadrant and 
started using the three-quadrant algorithm. 
Approximately 4 minutes after the beginning of the 
interference, the ORS voltage rapidly increases, quickly 
reaching its allowed maximum value of 4.0771 volts. It is 
saturated for about half a minute, then rapidly decreases, 
and stays at zero until the end of the interference. Note 
that pointing errors reach their maximum soon after the 
ORS voltage is reset to zero and are practically back to 
their usual range by the end of the interference period. A 
more detailed analysis is needed to determine the cause of 
this behavior of the ORS voltage. 


Finally, let us briefly analyze one of the numerous real- 
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EFFECT OF ROLL POINTING ERRORS 
ON TRMM YAW CONTROL 

The TRMM ACS 5 determines the yaw error by 
propagating the yaw angle between two DSS yaw updates 
with the gyro Z-rate. DSS measurements are used only 
twice per orbit when the Sun vector lies in the x B y B plane 
(in the FOV of either the DSS1 or the DSS2). This 
control law is valid 12 only under the assumption that there 
are no Earth pointing errors. To be more precise one 
needs to compute the Z-component of the angular velocity 

vector CDr in the geodetic coordinate system. The latter 
vector can be evaluated via the equation: 

“R = Abr m b -[0’“p-0] T . ( 14 > 

where A BR is the spacecraft attitude relative to the 
geodetic coordinate system, COg are gyro rates, and 
(Op is the pitch rate (±1 RPO). The yaw rate 
<|) y ~®R,z ’ correcte d for the roll pointing error, 
f , can thus be approximated as 


<t>Y “ 0) B,z + r03 P * (15) 

The integral over the first term in the right-hand side of 
this equation is commanded to zero, so that the yaw 

pointing error (Jjy can be represented as 
t 

<i> Y (t) = G>pj r(t')dt\ (16) 

0 

(It is assumed that misalignments of both DSSs have been 
accounted for and, consequently, that there is no yaw 
error immediately after the yaw update.) 

Fig. 7 compares the estimated yaw errors with their values 
propagated according to Eq. (16). One can see that 
propagated yaw errors are in excellent agreement with the 
estimation results. Thus, it is shown that any roll errors 
accumulated during ESA processing will contribute 
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Figure 7. Yaw Errors between DSS1 and DSS2 Yaw Updates 


significantly to large drifts in the yaw angle since roll- 
yaw coupling is not taken into account in the onboard 
propagation equation. 

CONCLUSIONS 

The present study reveals a systematic increase in 
penetration angle residuals during some parts of the 
TRMM orbit, which results in noticeable pitch and roll 
pointing errors. It is shown that the errors correlate with 
anomalous changes in the ORS voltages. This 
observation points to the ORS as a possible cause for the 
relatively poor pointing in this region. One advantage of 
this method is that an analysis of penetration angle 
residuals would make it possible to identify problems 
with the ORS for missions not having ORS voltages in 
telemetry (such as the series of NOAA satellites). 
Additional study is necessary to confirm this conclusion 
as well as to find the source of the observed anomalous 
behavior of the ORS heater. This study may also be 
important for the NOAA series, for which noticeable 
image distortions have been reported. 15 

Outside the regions of anomalously large pitch and roll 
errors, it is found that switching from the four-quadrant 


algorithm to the three-quadrant algorithm during periods 
of Sun or Moon interference in one of the ESA quadrants 
should keep Earth pointing to within the attitude accuracy 
specification. It is also demonstrated that relatively large 
roll and pitch pointing errors during interference periods 
should be attributed to the ORS voltage phenomenon 
discussed above (rather than to a lower accuracy of the 
three-quadrant algorithm). 

Direct computation of roll-yaw couplings prove that a 
large error in roll pointing plays a dominant role in the 
observed large drifts of the yaw angle between yaw 
updates using DSS measurements. This implies that one 
can use the spherical Earth approximation, with the only 
difference being that the infrared horizon, corrected for 
Earth oblateness, becomes latitude-dependent. 11 
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